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Abstract 

We investigate the effects of accretion of phantom energy onto primordial black holes. 
Since Hawking radiation and phantom energy accretion contribute to a decrease of the mass 
of the black hole, the primordial black hole that would be expected to decay now due to the 
Hawking process would decay earlier due to the inclusion of the phantom energy. Equiva- 
lently, to have the primordial black hole decay now it would have to be more massive initially. 
We find that the effect of the phantom energy is substantial and the black holes decaying now 
would be much more massive — over 10 orders of magnitude! This effect will be relevant 
for determining the time of production and hence the number of evaporating black holes 
expected in a universe accelerating due to phantom energy. 
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1 Introduction 



Numerous astrophysical observations are consistent with the standard cold dark matter 
model with the inclusion of an effective cosmological constant. A classical cosmological 
constant is generally avoided as quantum gravity attempts lead to a natural expectation 
of a cosmological term of Planck scale, which is totally at odds with the value required 
by observation (see for example [1]). As such, it is generally assumed that there is some 
physical field that comes into play after the Planck era and there is some principle that 
excludes the induced Planck energy cosmological term (see for example p]). This exotic 
field is often called "dark energy". According to the generally accepted modeling, the 
latter constitutes more than 70% of the total energy density of the universe while the 
matter component carries most of the remaining part [31 HI El El [7] . It is not clear how 
seriously to take the "quantum gravity" requirements considering that there is no viable 
theory of quantum gravity to date [UEj- Further, there is no clear evidence that the onset 
of quantum gravity will be at Planck scale and not orders of magnitude away from it [H] . 
There are alternate (classical) explanations of the observations available in the literature 
(for example by Wiltshire [TU]) but here we shall follow the generally accepted view of 
some form of dark energy providing the observed acceleration of the Universe. 

The observable universe locally appears to be spatially flat with an equation of state 
(EoS) parameter (the ratio of pressure to the energy density) uj{= Px/Px) — —1- The dark 
energy is an exotic vacuum energy with negative pressure and positive energy density 
which arises due to quantum vacuum fluctuations in spacetime. Caldwell and co-workers 
pn in\ considered the possibility of dark energy with the super-negative EoS parame- 
ter a; < — 1, which they called 'phantom energy'. It also gives negative pressure. The 
motivation to consider phantom energy as the candidate for dark energy arises from the 
observational data of the cosmic microwave background power spectrum and supernovae 
of type la. The phantom energy violates the general relativistic energy conditions in- 
cluding the null and dominant ones |T3] . Its implications in cosmology give rise to exotic 
phenomena like an imaginary value of the sound speed, negative temperature, the di- 
vergence of the scale factor a{t) and the energy density px ~ a"^^^"*"^-*, at a finite time 
resulting in a 'big rip', an epoch when the spacetime is torn apart [HI [151 US [13, HE] 
(see [19] for a review on the big rip singularity). However there are some attempts made 
recently in which the occurrence of a big rip singularity is avoided by phantom energy 
decay into matter [201 Ell [22] • Another attempt is the 'big trip', in which a wormhole 
accretes phantom energy and grows so large that it engulfs the whole universe [22] • A sim- 
ilar scenario is also proposed for black holes whereby the black hole event horizon inflates 
to swallow up the cosmological horizon, resulting in a naked singularity Moreover, 
quantum gravitational effects (if they exist) may avoid the big rip singularity. If the big 
rip cannot be avoided, the smaller the parameter uj, the closer the big rip will be to the 
present time. 
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Another weird property of phantom energy is that its accretion onto gravitationally 
bound structures resuhs in their dissociation and disintegration in a rather slow process. 
It was first analyzed in [25] for several gravitational systems like the solar system and the 
Milky Way galaxy. Initially, this possibility was investigated for a Schwarzschild black 
hole by Babichev et al [2S] , who showed that the black hole mass goes to zero near the big 
rip. Interestingly, in this scenario larger black holes lose mass more rapidly than smaller 
ones. Later on, their model was extended to the Reissner- Nordstrom [28], Schwarzschild 
de Sitter [29] and Kerr-Newmann [301 EI] black holes. It should be mentioned that it has 
been argued that the mechanism of accretion followed by Babichev et al is stationary and 
does not possess the shift symmetry [27] and hence that the mechanism of dark energy 
accretion is not realistic and consistent. Nevertheless, we shall follow the Babichev et al 
analysis, leaving the detailed analysis for subsequent work, as the effect will be technically 
difficult to apply and we believe will not make a substantial difference for phantom energy 
in the neighbourhood of a primordial black hole. It has also been argued that when the 
back-reaction effects of the accretion process are included in the analysis of Babichev et al 
[25] , the black hole mass may increase instead of decreasing [321 E3] , thus avoiding the big 
rip. Also in cyclic cosmological models, black holes do not tear apart near the turnaround 
but preserve some nonzero mass [Ml IH3] . We shall ignore the big rip issue here. 

2 Hawking evaporation of black holes 

We are interested in studying the effects of accretion of phantom energy on a static 
primordial black hole. Carr and Hawking [3S] in 1974 considered the formation of black 
holes of mass lO^kg and upwards in the early evolution of the universe. After their 
attempt, several authors investigated various scenarios of PBH formation EHl EHl SO]- 
The existence of these small mass black holes was based on the assumption that the 
early universe was not entirely spatially smooth but there were density fluctuations or 
inhomogeneities in the primordial plasma which gravitationally collapsed to form these 
black holes. Unlike the conventional black holes that are formed by the gravitational 
collapse of stars or mergers of neutron stars, the primordial black holes (PBH) are formed 
due to the gravitational collapse of matter without forming any initial stellar object. The 
mass of a PBH can be of the order of the particle horizon mass at the time of its formation 



Therefore PBHs that formed in the early history of the universe must be less massive 
while those that formed later must be more massive. Black holes formed at Planck time 
10~^'^s would have Planck mass 10~®kg. 

Using classical arguments, Penrose and Floyd showed that one can extract rotational 
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energy from a rotating black hole ^3]. Penrose went on to argue (see [H] and references 
therein) that one could take thermal energy from the environs of a black hole and throw 
it into the black hole to get usable energy out. This would apparently reduce the entropy 
around the black hole. As such, he had argued that there must be an entropy of the 
black hole that increases at least as much as that of its environs decreases. Hawking had 
pointed out that in any physical process the area of a black hole always increases [H] 
just as entropy always increases. This led Bekenstein [15] to propose a linear relationship 
between the area and entropy of a black hole. Thus Bekenstein ^61 W7\ generalized the 
second law of thermodynamics to state that the sum of the entropy of the black hole and 
its environs never decreases. However, at this stage it seemed that the connection between 
black holes and thermodynamics was purely formal. At this stage Fulling pointed out that 
quantization of scalar fields in accelerated frames gives an ambiguous result [48j, which 
seemed to yield radiation seen in the accelerated with a fractional number of particles. 
Hawking repeated the calculation for an observer near a black hole and obtained the same 
result by various methods and found that the radiation had a thermal spectrum |39j . This 
led him to propose that mini-PBHs would evaporate away in a finite time 



The corresponding Hawking evaporation process reduces the mass of the black hole by 



dM 



dt 



(2) 



where a is the spin parameter of the emitting particles. Integration of Eq. (2) gives the 
evolution of PBH mass as 

/ tV^ 

Mur = MAl- — \ , (3) 
where the Hawking evaporation time scale thr is 

It is obvious from Eq. (3) that as t — > thr, the mass Mhr — ^ 0. Plugging in thr = to (the 
current age of the universe) in Eq. (4) gives the mass lO^^kg of the PBH that should 
have been evaporating now. Hence from Eq. (1), it can be estimated that these PBHs 
were formed before about 10~ sec. For Mi > lO^^kg, a = 2.011 x 10"^, hence Eq. (4) 

implies thr ~ 2.16 x lO^^* (^g) sec. While for 5 x lO^^A;^ < M, < 10^\g, a = 3.6 x 10"^ 

then Eq. (4) gives thr — 4.8 x 10^^*^ (^) therefore detecting PBHs would be a 

good tool to probe the very early universe (closer to the Planck time). The evaporation 
of PBHs could still have interesting cosmological implications: they might generate the 
microwave background ^52j or modify the standard cosmological nucleosynthesis scenario 
[53j or contribute to the cosmic baryon asymmetry [S^ . Some authors have also considered 
the possibihty of the accretion of matter and dust onto the seed PBH resulting in the 
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formation of super-massive black holes which reside in the centers of giant spiral and 
elliptical galaxies [55] . 



3 Phantom energy accretion onto black hole 



The FRW equations governing the dynamics of our gravitational system are given by 

H^={z] = —iPm + Px), (5) 

I 

-[Pm + Px(l • (6) 



a J 3 

d AttG, 
a 



Here pm and denote the energy densities of matter and the exotic energy densities 
respectively. The scale factor a(t) goes like [17] 

a{t) = ^ {t>to), (7) 

[-u; + (1 + uj)t/to\ 3(1+-) 

where to is the time when the universe transits from matter to exotic energy domination 
(which is roughly equal to the age of the universe). Notice that the scale factor a{t) 
diverges when the quantity in the square brackets in Eq.(7) vanishes identically i.e. 



t* = T^to . (8) 

1 + UJ 



Subtracting to from Eq. (8), we get 



t*-to = -^to . (9) 

I + 00 

The evolution of energy density of the exotic energy is given by 

p-'=6nGil + ujrit*-ty. (10) 
A black hole accreting only the exotic energy has the following rate of change in mass 



dM{t) 



dt 



2 



M\p^+p,). (11) 



It is clear that when Px + Px < 0, the mass of the black hole will decrease. We are 
particularly interested in the evolution of black holes about and after t = to since the 
dark energy is presumably negligible before that time and may not have any noticeable 
effects on the black hole. Using Eqs. (9) and (10) in (11), we get 



dM{t) 



dt 



8GM^ , 
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Therefore the mass change rate for a black hole accreting pure exotic energy is determined 
by Eq. (12). For the phantom energy accretion, the time scale is obtained by integrating 
Eq. (12) to get 



Mit) = 1 - 



t 

where tx is the characteristic accretion time scale given by 

.1 _ IGttG^ 



c5 



-Mi{px+Px) ■ 



Using Eqs. (9) and (10) in (14), we get 



3c3 tl 
8GMi{l + u;) 



(13) 



(14) 



(15) 



4 Evolution of mass due to phantom energy accretion and Hawk- 
ing evaporation 



The expression determining the cumulative evolution of the black hole is obtained by 
adding Eqs. (2) and (12) i.e. 



dM{t) 



dt 



dM 



Total 



dt 



+ 



dM 



hr 



dt 



We write the above equation as 



where 



M2 3c3 tl 



— = -aM^ - — , 
dt M2' 



8G e hc^a 



Here e = —u — 1. Thus (18) can be written in the form 

If M^dM 

- / dt 



b J 1 + f 



To integrate above equation, we assume 



X 



1/4 



(16) 
(17) 

(18) 
(19) 



(20) 
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which yields 
We note that [56 



- / dt 



x^dx 



{a%Yl^ J l + x 



(21) 



l + x 



2n 



2n 



E 

k=l 



COS 



2n 



In 



2xcos 



2A; - 1 
2n 



■K + X^ 



1 /m7r(2k — li \ ^ 
- 2^ sm f ^ ] tan"^ 



n 



k=l 



2n 



X - cos I ^ ITT 



Sm I I TT 



2n 



m < 2n. (22) 



In our case, m = 3 and n = 2, hence the above equation yields 



x'^dx 



In 



1 - V2x + 



1 + V2x + x2 



1 + 4v^ 
On substituting the value of x above, we obtain 



1 _i / V2x 
tan 



2^/2 



l + x' 



(23) 



t = U 



1 



4^2 



In 



1-V2(t 



1/4 



M+ (f 



1/2 



/ \ 1/4 / \ 1/2 



+ 



2^2 



tan 



-1 



1/4 

2(f) M 



l+lf 



1/2 



M2 



(24) 



We now redefine the values of a and h by assuming M = mMi, where m is a dimensionless 
parameter and Mj is the initial mass of the black hole. Thus (18) becomes 



dm 
~dt 



I 2 

-a m 



h' 



(25) 



where a' = aMi and h' = h/Mf. For the terms to be equal strength, we require a' ~ h' . 
Thus 

M,^(-] . (26) 



Now 



We can normalize 



b _ 3hcHla 
a ~ 8G^e ' 



1/4 



?T 



a/ 



or, e 



1/2 



3nc^t? 



(27) 



V2{ % 



1/4 



Replacing p = a' /h' 



Mf 



Mf (the ratio of the phantom component to the 



Hawking component, in the energy radiated) the above equation becomes 



t = tn 



1 - 



l-V2pl/4+pl/2 



l + V^pl/4+pl/2 



+ 2v^ [ l+pi/2 



(29) 



Moveover, the power emission due to Hawking evaporation from the stationary black hole 
of mass M > lO^^g 



P = 3.458 X 10^''(M/g)-2ergs-\ 
and for mass 5 x lO^^g ^ M ^ 10^''g, 

P ^ 3.6 X 10^^(M/10^^g)"2ergs-^ 
In our analysis, the mass in the above two expressions is replaced by 



M 



( 8G^e 



(30) 



(31) 



(32) 



Now choosing e = 0.1, we obtain M = 8.74029 x lO^^g which will be evaporating now due 
to the combined effects of phantom energy and Hawking radiation. Then using Eq. (30), 
the corresponding power emission will be P = 4.52661 erg s~^. We can compare this result 
with that of a black hole of mass M ~ 1.05 x lO^^g evaporating just now due to Hawking 
radiation only. The corresponding power emission will be P ~ 3.144 x lO^^erg s~^. Note 
that the power emission from a black hole decreases when the effects of phantom energy are 
incorporated. Similarly, for very large values of e ~ 10^^ would give M = 2.763923 x lO^^g. 
Using this mass in (31), the power emitted is 4.52661 x 10^'^erg s~^. However, such large 
values would lead to a very early big rip and hence must be excluded. Thus black holes 
~ lO^^g are of more interest for observational purposes since these are the ones that should 
be evaporating now. 



5 Conclusion 

In this paper, we have analyzed the Hawking radiation effects combined with the phantom 
energy accretion on a stationary black hole. The former process has been thoroughly 
investigated in the literature. However there is as yet no observational support to it. 
According to standard theory it is assumed that after the formation of PBHs (of mass 
~ lO^^kg with a Hawking temperature lO^^K), they would absorb virtually no radiation 
or matter whatsoever during their evolution and radiate continuously till they evaporate 
in a burst of gamma rays at the present time. This scenario assumes that the Hawking 
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temperature for such black holes was always larger than the background temperature of 
the CMB. Strictly speaking, this cannot be true. Consequently PBHs could have accreted 
the background radiation (and even some matter) and grown in mass. Hence there should 
be no PBH left to be evaporating right now However, the above scenario is modified 
when phantom energy comes into play. When phantom energy and the Hawking process 
are relevant the total life time scale of the PBH is significantly shortened and the formation 
of the PBH exploding now is delayed. 

From Eq. (29) we obtain the time as a function of mass instead of getting mass as 
a function of time. To make sense of the results we need to obtain the evolution with 
time. This is done by inverting the explicit function. We have plotted the normalized 
time r = t/to against the dimensionless mass parameter m and m against r for different 
choices of the parameter p, in Figures 1 - 10. It is observed that increasing p increases 
the steepness of the curve specifying the mass evolution. Therefore the black hole loses 
mass faster for larger p till it vanishes at r = 1, the present time. In particular. Figures 7 
and 9 show the same evolution of mass for larger values of p. It appears that the graphs 
contain a redundant (or nonphysical) part of the mass evolution and the only physically 
interesting section is above the horizontal curve crossing t = 0. Thus in effect, see Figures 
8 and 10, the initial mass of the black hole must be taken 0.45Mj of the value given by 
for p = 5 and about 0.315Mj for p = 10. It is obvious that the results are very insensitive 
to changes of the parameter e for the phantom energy. As such, they can be regarded as 
fairly robust. 
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Figure 1: The normalized time is plotted against the mass parameter for = 0.1. 
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Figure 2: The mass parameter is plotted against the normalized time for p = 0.1. 
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Figure 3: The normalized time is plotted against the mass parameter for p = 0.5. 
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Figure 4: The mass parameter is plotted against the normalized time for p = 0.5. 
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Figure 5: The normalized time is plotted against the mass parameter for p = 1 
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Figure 6: The mass parameter is plotted against the normalized time for p = 1 



15 



0.2 0.4 0.6 O.S 1 

111: 



Figure 7: The normalized time is plotted against the mass parameter for p = 5. 
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Figure 8: The mass parameter is plotted against the normalized time for p = 5. 



16 



Q ''d'^2-_ Q.4 -oVe' 0.8 .i 



Figure 9: The normalized time is plotted against the mass parameter ior p = 10. 
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Figure 10: The mass parameter is plotted against the normalized time for p = 10. 
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